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We calculate the magnetization of a two-dimensional superconductor in a perpendicular magnetic 
field near its Kosterlitz-Thouless transition and at lower temperatures. We find that the critical 
behavior is more complex than assumed in the literature and that, in particular, the critical mag- 
netization is not field independent as naive scaling predicts. In the low temperature phase we find a 
substantial fluctuation renormalization of the mean-field result. We compare our analysis with the 
data on the cuprates. 



I. INTRODUCTION 

The study of fluctuation superconductivity received a 
tremendous impetus from the discovery of the cuprate, 
high temperature superconductors 1 -. It was quickly real- 
ized that these were materials with large Ginzburg pa- 
rameters with the additional feature that many of them 
are also highly two-dimensional. More recently, interest 
has focussed on the systematics of their superconducting 
properties with doping and it has been found that this 
anisotropy is further enhanced with undcrdoping 2 . Con- 
versely, the extent to which the anomalous properties of 
the underdoped compounds in the "pseudogap" region 
can be attributed to superconducting fluctuations is a 
question of considerable interes^*iiSiSiL2iS. 

In this context, we report here a study of the mag- 
netization of a two-dimensional superconductor, or more 
precisely of a stack of decoupled two dimensional super- 
conducting layers, in a perpendicular magnetic field H. 
By combining Kosterlitz-Thouless renormalization group 
flows and explicit computations for plasmas, we find the 
field and temperature dependence of the magnetization 
density, M(H,T), for temperatures T near to or below 
the Kosterlitz-Thouless transition temperature Tkt at 
fields H c i <C H <§C H C 2- These results are interesting 
on three immediate fronts. First, the rough magnitude 
of M(H, T) and the trends exhibited in Fig. 1 are in 
accord with much of the data on the strongly layered 
cuprates in the field regime where a single layer theory 
is expected to apply. For example, this holds for the re- 
cent data of Wang and collaborators 4 which have demon- 
strated the existence of substantial fluctuation magneti- 
zation in the same region of the phase diagram of BSCCO 
which was previously observed to exhibit a large Nernst 
signal 5 >. Second, we find that M is not independent of H 
at criticality, as assumed in earlier analysis of the "cross- 
ing point" phenomenon in studies of the cuprates 1 ^. In- 
stead, the critical magnetization follows 



M : 



log 7i log 



(1) 



where d is the distance between layers, </>o is the flux 
quantum, H is the magnetic field component perpendic- 
ular to the layers and do is the microscopic short-distance 
cutoff length (71 and 72 are constants). Third, we find 



that below Tkt, M exhibits a low field growth indicative 
of the expulsion of vortices from the KT phase, 
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where p s is the 2D superfluid density, related to the in- 
plane penetration depth via p s — d(pQ / {Am 1 A 2 /xq) (73 is a 
constant that vanishes as T approaches Tkt, see below). 
Note that this expression gives a simple relation between 
the superfluid density and the derivative of the magneti- 
zation with respect to the logarithm of the field, that can 
be used to determine when the system is indeed behav- 
ing as a stack of essentially decoupled two-dimensional 
Kosterlitz-Thouless films. 

These results and the limits of their applicability to 
actual layered superconductors are derived in Section II, 
with some details relegated to an Appendix. In Section 
III we consider experimental signatures of our analysis 
and the existing experimental situation. We close with a 
brief summary and some open questions. 



II. THEORY 

A. Definition of the problem 

The average magnetization density of a superconduc- 
tor in an external field H is 



At0 M = B(H) - /UoH, 



(3) 



where B(H) is the average (uniform) magnetic field, com- 
puted via 



df, 



(B) 



dB 



H, 



(4) 



where f scrn {B) is the total free energy density of super- 
conducting matter coupled to fluctuating magnetic fields 
with average field fixed at B. Let us define 
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In type II materials and for H H c i field energy is 
well-approximated by its uniform value, and therefore, 
following the equation above, the magnetization is to be 
computed from the free energy of a charged superfluid in 
a uniform external magnetic field. 

We consider a stack of decoupled two dimensional lay- 
ers. Hence we will replace the three dimensional free 
energy density / SC (B) by \f(B) where f(B) is the free 
energy density of a single layer in the presence of a mag- 
netic induction B, now restricted to be perpendicular to 
the layers, whose effect will be to impose a density B/(po 
of field- induced vortices. 



B. Coulomb gas formulation 

To compute / and thence M we shall resort to the 
standardii mapping of the two-dimensional vortex prob- 
lem onto a two component Coulomb plasma whose 
Hamiltonian is given by (see Appendix A) 



nr. . 

H = N T E c0 + irp s0 ^2 PiPj lo S + H B - 



(7) 



l<3 



The number of vortices of charge pi = ±1 is N±. The to- 
tal number of vortices is Nt = N + + N- and is allowed to 
fluctuate by the addition and removal of neutral vortex- 
antivortex pairs, but the net charge Q = N + — N- = 
L 2 B/c/)q is constrained by the field B. E c0 is the bare 
vortex core energy, ao — the bare short distance cutoff 
(e.g., vortex core radius), p S Q is the bare superfluid stiff- 
ness (inverse dielectric constant in the plasma language), 
L 2 is the area of the system. A uniform background 
density of charge to make the system neutral is neces- 
sary to insure a proper thermodynamic limit. The total 
Coulomb interaction potential between the vortices and 
this background and between the background and itself 
is in the constant Hb that depends on B and on the size 
and shape of the system, but not on the vortex configu- 
ration. 



C. Renormalization and matching 

Standard Kosterlitz renormalization group (RG) 
methodsAiiiSiiii can be generalized to the present, non- 
neutral situation. Our strategy will be to renormalizc 
until we can match on to a solution by an exact or ap- 
proximate method. For convenience, the bare superfluid 
stiffness and core energy are captured by the bare cou- 
plings x = 1 - jj^, and yo = 27re~ E <= o/ ' csT , which 
we will renormalize. In addition to these couplings we 
shall be interested in computing the free energy den- 
sity, so we keep track of the configuration-independent 
term C in the renormalized Hamiltonian, generated as 
degrees of freedom are integrated out, i.e., at any inter- 
mediate step in the RG process the partition function is 
Z = er L c / a Tre _/3ff , where a is the renormalized cutoff. 



Note C is defined to be dimensionless. We also define the 
dimensionless density of field-induced vortices (number 
per renormalized cutoff area) as n — Ba 2 /4>a; its bare 
value is uq = Bo^/^q. 

Provided the total number of vortices and antivortices 
per cutoff area is small compared to one, the following 
differential equations describe the renormalization upon 
increasing the cutoff to a = a^b and integrating out neu- 
tral vortex-antivortex pairs with spacing less than a: 
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The last line shows the trivial renormalization of the 
number of field induced vortices per cutoff area. A some- 
what unexpected fact about these equations is that the 
presence of field induced vortices leaves the zero- field flow 
equations (Eqs. ijHJ), I© and (|10Jl ) unaffected (also see Sec. 
IV). As we shall see, this simplifies the calculation of the 
magnetization. 

By straightforward integration one obtains, with 
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C(b) = b 2 C(0)-b 
n(b) = n^b 2 
with coshO< : 



c cotanh(9< + 2c log b) 
c cosech(9< + 2c log b) 
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^ for T < T KT , and the solutions, 



x>(b) 
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with cos9> = — for T > T KT . 

The free energy density of the original problem is then 
computed by running the RG to scale 6_r, where it can 
be expressed as a sum of two parts 



/ 



T R + k B TC{b R ) 



(14) 



where is the free energy per renormalized cutoff area 
of the vortices that have not yet been integrated out. 
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[Note: since we have chosen to rescale lengths in our 
RG, the free energy density / is computed by dividing 
the residual free energy, Tr, and C by the rescaled 
area, a\] The observation that field induced vortices 
do not affect the other flows, in particular that C has 
no explicit dependence on B, simplifies the calculation of 
the magnetization density: 
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thus reducing the problem to computing Tr. 

The renormalized density of field-induced vortices, n, 
grows under the RG flow. For T < Tkt thermally in- 
duced vortices become increasingly dilute and the sys- 
tem approaches a one component plasma in the low- 
field limit where the flow can continue to large b^. For 
T > Tkt, the renormalized density of thermal vortices 
first decreases and then increases, so that the system is 
always a two-component plasma with a scale-dependent 
charge imbalance. Our principal results for the asymp- 
totic behavior of magnetization (Eqs. (1) and (2)) follow 
from the exactly known free energy of a dilute one com- 
ponent plasma^* whose parameters, the charge density, 
core energy and dielectric constant, are given by the RG 
flows above. A Debye-Hiickel approximation for the free 
energy of the two-component plasma reproduces the cor- 
rect low field behavior for T < Tkt and provides a rea- 
sonable approximation elsewhere, including T > Tkt', we 
shall use it to generate plots of the magnetization versus 
field in the vicinity of the transition. 



D. Results 

1. Low fields, T < T KT 

For T < Tkt the above Kosterlitz-Thouless RG flows 
go to small fugacity y, where the leading terms in the 
exact free energy are known for low density (small field). 
Thus we run the RG to a matching scale an = bRCio, 
where the parameters are ur — Ba 2 R /4>Q, xr and dr. 
We stop and match when n, which is growing, becomes 
of order y. The free energy density of the remaining 
vortices, to leading order in the small parameters ur and 

■ 14 

VR 

J 7 R = n R k B T (log— +x R \ognR + 0(l)) . (16) 

VR 

The first term is the (renormalized) core energy of the 
field-induced vortices, while the second term contains 
their entropy and (renormalized) interaction energy. The 
thermally-excited vortex-antivortex pairs do not enter at 
this order. Thus, from (| 1 51) the magnetization is 

M = -^-{\og—+XR\ogn R + 0(l)). (17) 
d<po VR 



To start, we consider the regime of small fields below 
Tkt- For b\ ^> e 2j4 ^ c , where 74 is a constant, the renor- 
malized parameters are computed from the Kosterlitz RG 
equations as (cf. Eq. I|12|) ') 
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so the low field magnetization is 
M 



ksT <f> 



Bag 



log ^ ±M + 
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in agreement with our Eq. (2) above (73 ~ c 1 ^). The low 
field regime of validity of this expression requires ur <C 1 , 
which, with the above constraint on 6r, translates to 



B « %- 2 ^ 



(21) 



Observe that this condition defines a crossover length 
scale 



a e 



74/c 



(22) 



which has the functional form of the correlation length 
above the transition, and yet is defined — like the Joseph- 
son correlation length in Goldstone phases — below the 
transition. 

For higher fields and temperatures near Tkt, there is 
a crossover to the critical behavior (1), that we derive 
below. Note that the slope of M vs. logB at low field 
vanishes linearly in c as T approaches Tkt from below, 
while the crossover field below which this is the behavior 
vanishes exponentially. 

If, on the other hand, one is near the Kosterlitz- 
Thouless transition (so c ^ 1), the RG flows for 1 <C 
b% < e 2 ^' c follow 



(23) 
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Here when we match at n R ~ yR, it is at yR, 



log(Ba%/4> ) 



and the resulting magnetization is 
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which is the same as Eq. (1). The field range where this 
result applies is 



(I; 



2 > B > ^2" 



(25) 



< 



which is an intermediate field range for T near to but be- 
low Tkt, and is all low fields for T = Tkt- Our method 
of analysis in principle also gives the functional form of 
the crossover scaling function between these two regimes, 
but we have not been able to write this function in any 
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concise form. However, below we obtain an approximate 
scaling function covering both regimes. 

It is worth noting that in the vicinity of Tkt the 
matching, which we choose to do at ur = y Rl occurs 
at \xr\ -C 1, so it is the renormalized core energy (the 
first term in Eq. I|17|l) that is dominant in determining 
the magnetization. 



2. Approximate results, 



Above Tkt the RG flows at asymptotically small fields 
are terminated at a finite fugacity and a$b ~ £ (the zero 
field correlation length), so that as B — ► the magneti- 
zation vanishes as 
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This expression is consistent with an earlier linear re- 
sponse result^. 



Building on this last observation we now consider a 
Debye-Hiickel mean field theory whereby the matching 
free energy is approximated by the contributions from 
the renormalized core energy and the entropy. This will 
allow us to obtain numerical results for the magnetiza- 
tion which are consistent with our previous exact results 
and can be extended to above Tkt- However, we must 
caution that these are approximate results, and the pre- 
cise quantitative values of the magnetization do depend 
on aspects of the calculation that are not a priori con- 
strained by what we already know about these materials. 
These unconstrained freedoms in the approximation we 
present below include (i) the form of the matching free 
energy, in particular one might also include an approx- 
imation to the energy of vortex-vortex interactions, (ii) 
the scale at which the matching is done, and (iii) the 
bare vortex core energy. 

Specifically, we approximate the residual free energy 
density at matching as 

•F R = {n% + n B )E cR + k B T{n^ logn^ + n R logn^) (26) 

where vortices and antivortices remaining at matching 
are at densities and n R , respectively. The field 
constrains the difference between these ^densities to be 
lift — rift = tir, but their mean, Wr = , can vary 

and will take on the value that minimizes the free energy: 

O 'T" 2 

— ^ = 2E cR + k B T\og[n^ 2 - ^R] + 2k B T = . (27) 
(Mr 4 

The matching condition is defined by asking that the 
average density of vortices be 1: 
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and this is used to solve for or{B,T). The choice of 
precisely unit density is clearly arbitrary — we make it for 
concreteness. With these assumptions the magnetization 
is given by 



3. Comments 
Finally, three comments are in order. 

First, the divergence of M below Tkt signals the ex- 
pulsion of vortices from that phase — indeed the logarith- 
mic divergence of the free energy of a single vortex with 
system size. The divergence we find at Tkt likewise sig- 
nals the now weaker expulsion of vorticity; At Tkt the 
free energy of a single vortex grows with system size as 
log(logi). 

Second, the divergence of M(B) at low B is a correct 
statement about the magnetic induction dependence of 
the free energy of a single layer in the limit of infinite pen- 
etration depth. It does not imply the divergence of an ac- 
tual three dimensional magnetization density. For a stack 
of planes, the self-consistency implicit in Eq. (6) keeps 
M(fj,oH) from blowing up — instead we get the Meissner 
phase (B = 0) at sufficiently low applied fields. Our 
results are valid in the experimentally relevant regime 
M <C H (equivalently H ^S> H c x) where this distinction 
is not important. 

Finally, for T^Tkt as we have already remarked B* « 
separates the small field behavior from the critical 

log log 1/B dependence at larger fields. Interestingly, the 
magnetization displays a noticeable temperature depen- 
dence even at these larger fields. It is noteworthy that 
this dependence is opposite to that implied by the ex- 
plicit prefactor of temperature (e.g., in Eq. (J20J)), which 
by itself predicts an unintuitive enhanced diamagnetism 
at higher temperatures. To understand this important 
detail consider raising T in the vicinity of Tkt- This 
leads to an increase in vr, which in turn gets strongly 
amplified by the singularity at vr = of the functional 
dependence M ~ logl/ 'y R . This variation always over- 
whelms the contribution from the prefactor of k B T,so 
the combined effect is a reduction of \M\ upon raising 
the temperature, as expected. 
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which is plotted in Figs. ^ 121 ancl CH 



(29) E. Dimensional crossovers in layered and quantum 

systems 

(30) 

In sufficiently anisotropic materials we expect that 
various weak interlayer couplings present will only af- 
fect the physics at long length scales and can therefore 
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FIG. 1: Field dependence of magnetization for a range of 
temperatures between T = 77 K and T = 83K, with the 
dashed curve corresponding to T = Tirr = 80A" (others are 
spaced in 0.5if increments). Here we take E c = -Kp 3 o(T), 
which then implies ksTuT ~ 1.13p»o (Tkt)- We chose 
a = 30^, d = 15A, so that Bcfi/fa « 0.04 at 10 Tesla and 
IzbTkt / {d<j>o) m S50A/m. Finally, we model the experimen- 
tally observed temperature dependence of the bare superfluid 
stiffness by p s0 (T) / p s0 (0) = 1 - T/12QK (see, e.g., Ref.©. 
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FIG. 2: Semi-logarithmic plot of magnetization. Parameters 
are as in Fig. 1. 



be treated to simply restrict the domain of applicabil- 
ity of the purely two-dimensional treatment above. As 
the application of the magnetic field introduces a new 
length scale, the intervortex separation or the magnetic 
length, this restriction means our two-dimensional be- 
havior will cross over to three-dimensional behavior near 
a low crossover field. There are two basic types of inter- 
layer interactions here: Josephson couplings induced by 
charge fluctuations transverse to the layers, and electro- 
magnetic couplings between currents in different layers 
mediated by the fluctuating magnetic field. We consider 
their effects in turn. 

The Josephson coupling, J, will a) induce true three- 



dimensional long range order in the low temperature 
phase in zero field and alter the universality class of the 
transition to that of the 3DXY model, b) shift the transi- 
tion temperature upwards to T c and c) cause the vortices 
to crystallize three dimensionally at any temperature be- 
low T c for fields less than the melting field B m (T). While 
this is a problem that has been studied in some detailii 
we can get a feeling for the scales involved by considering 
renormalization of Ja 2 , the Josephson coupling per cut- 
off area, near the (decoupled layer) Kosterlitz-Thouless 
fixed point, 



d(Ja 2 ) 
d log b 



(2 



)(Ja 2 ). 



(32) 



From this equation we can obtain the scaling of both the 
shift of the critical temperature in zero field, ST C , and 
the field scale, Bj, marking the crossover between two- 
and three-dimensional behaviors. This is done by eval- 
uating the length scale at which the renormalized value 
of Josephson coupling is comparable to the in-plane stiff- 
ness, Ja 2 w p s (which we shall set to its critical value 
p s = DzbTkt I ~k , since its renormalization is compar- 
atively less important than that of J a 2 ). The upward 
shift of the critical temperature due to the Josephson 
coupling, 8T Cl can be estimated by identifying this scale 
with the zero-field correlation length, giving 
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where 75 is a non-universal constant related to the be- 
havior of c near Tkt via 75c « yj (T — Tkt) /Tkt, 76 is 
another constant of order one, and Jo is the bare Joseph- 
son coupling per unit area. Similarly, by substituting the 
magnetic length in place of £, we find for the crossover 
field at T KT , 



Bj = 



>oa 



kBTKT J 



(34) 



The case of electromagnetic coupling alone is a lit- 
tle muddier. The situation in finite magnetic fields is 
that the vortices in different layers now experience an 
attraction and thus can crystallize three dimensionally 
even when a single layer is a vortex liquid. Naively, the 
crossover field is of order the melting field, which is esti- 
mated to bei*i& 



i>o ( \i> yk B TKTa\ 
<j> 2 d 



(35) 



where A is the penetration length at Tkt- The impact 
of these magnetic interlayer couplings on the zero field 
transition does not appear to be a settled problem. The 
interlayer interactions are also logarithmic and thus are 
marginal operators at face value. This has led to asser- 
tions that the actual transition is still in the Kosterlitz- 
Thouless universality classi However, a renormalization 
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group analysis by Timrr. 19 finds flows at variance with 
this scenario. For our present purposes, we note that this 
crossover field is, for the case of large k that we are inter- 
ested in, less than /iqH c i] it is this crossover that gener- 
ally restricts the validity of our two-dimensional analysis 
toB» fioHcx. 

The scales Bj and Bb mark the rough boundary be- 
tween two and three dimensional physics. At Tkt the 
low-field three dimensional state is crystalline, and the 
high-field two dimensional dependence of M crosses over 
to fairly standard (but with renormalizcd parameters) 
low-field behavior in the Abrikosov vortex lattice. At 
higher temperature at the three dimensional T c , in the 
presence of Josephson couplings, the crossover is instead 
to the M ~ —\[iL behavior expected at low field in the 
3DXY critical regime. 

Finally, one additional crossover, this time at large 
fields, is possible when the thermal phase transitions take 
place in proximity to a quantum phase transition out of 
the superconducting state — as may be germane in the 
case of the underdoped cuprates. Standard scaling, when 
applied to magnetization of the 2 + 1-dimensional quan- 
tum critical theory predicts M ~ —y/~H, with crossover 
to this behavior taking place at sufficiently large fields 
where quantum fluctuations of the order parameter are 
important. The transition between the forms derived in 
this paper and this regime would be a striking signature 
of such fluctuations. 



strong amplitude fluctuations at high fields near mean- 
field H c2 (see, e.g., Ref. |21|). Calculations of this kind 
have no overlapping regime of validity with ones pre- 
sented above for the low field regime. It would be worth- 
while to see if the two sets of results can together capture 
the behavior over the entire range of fields of interest. 

Finally, Bulaevskii, et <d~ have considered the ef- 
fects of thermal fluctuations on the magnetization of an 
Abrikosov lattice in a layered superconductor with the 
Josephson coupling being the dominant source of the 
three-dimensional order. Treating the phonons of the lat- 
tice to quadratic order they compute the entropic correc- 
tion to the Ginzburg-Landau- Abrikosov free energy. At 
low fields B <C B cr they find a field independent correc- 
tion to the leading (Abrikosov) logarithm. For B 3> B cr 
they report a correction which is itself a logarithm and 
leads to an expression identical to our Eq. (2). The scale 
B cr is the Lindemann estimate for the melting field of the 
vortex lattice they begin with, so the computation is not 
really valid in the regime of the logarithmic dependence 
where the system is now a vortex liquid. Nevertheless, 
at sufficiently low fields B <C 0o/£< and T < Tkt their 
result and ours agree, a sign that the vortex liquid is lo- 
cally quite similar to the crystal for these parameters. At 
higher fields and at Tkt, however, their reasoning breaks 
down as our calculation explicitly demonstrates. 

III. EXPERIMENTS 



F. Prior work 

Early experiments on the most anisotropic cuprates 
saw a near crossing point for the curves M(T) taken at 
varying fields B. This is equivalent to the statement that 
they found a temperature at which the magnetization 
was field independent. This was interpreted in the scal- 
ing framework 2 as evidence for two dimensional critical 
behavior with the one-parameter scaling form, 



f(t,B) 



\hBe) 



(36) 



for the free energy density. As we have shown in this pa- 
per, the scaling is not so simple at the Kosterlitz-Thouless 
transition, with its two marginal operators, and magneti- 
zation instead has a weak double-logarithmic dependence 
on the field. Aside from this general scaling argument, 
there are three types of prior computations that we are 
aware of (not listed in chronological order): 

First, Gaussian fluctuations for the Ginzburg-Landau 
theory in d = 2, yield a field independent magnetization 
at the Gaussian (mean field) transition temperature^. 
This is consistent with the absence of any marginal opera- 
tors at this unstable fixed point. Our calculation replaces 
this result as an account of the true low-field critical scal- 
ing of M, 

Second, the so called "lowest Landau level" approxi- 
mation has been used by Tesanovic and others to study 



We now turn to the existing data on the cuprates, the 
systems that have motivated this work, and some sug- 
gestions for experimental tests of the theory. 

Many of the cuprates are highly two dimensional and 
appear to become increasingly so with underdoping. The 
two most commonly used diagnostics of anisotropy are 
the resistivity and superfluid density tensors. Using ei- 
ther of these in materials such as BSSCO-2212 we arrive 
at estimates of the anisotropy of order 10 -4 to 10~ 6 be- 
tween the ab (Cu-O) planes and the c-axis. A somewhat 
more direct measure of the anisotropy can be obtained 
from observations of the c-axis plasma resonance due to 
interlayer Josephson coupling. From the results of Ref. 
I23I for Josephson coupling per area Jo ~ 10 _8 Joule/m 2 
we obtain 
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(37) 



For such anisotropies and other parameters appropriate 
to the cuprates, e.g., as in caption of Fig. 1, we find 
8T C /T KT « 0.02, Bj fa 0.0002 Tesla and B B « 0.004 
Tesla using the estimates derived in Section HE. This 
indicates that our two dimensional theory should give a 
useful account of the superconducting fluctuations for a 
reasonable range of fields. 

The scale for magnetization effects is set by fc -^"^ T 
which is « 350^4/to for our parameters. The dimension- 
less factors multiplying it in our expressions are not too 



different from unity. This order of magnitude estimate 
is consistent with experimental observations 424-25 . The 
trends in the temperature and field dependence of the 
magnetization we show in our figures here are mostly 
in good qualitative agreement with those seen in the 
experiments^^ 5 -. The most detailed published inves- 
tigation is that by Kogan and collaborators 24 (sec also 
Ref. |25|) who were inspired by the predictions of Ref. I22I 
They reported evidence that M ~ log if with a coeffi- 
cient that changes sign at Tkt- While this claim below 
Txt is consistent with our analysis, this is not so above 
Tkt where the functional form is no longer a logarithm. 
This suggests that a reexamination of that regime is in 
order. 

Finally, we record the salient results of our analysis 
that can be tested against careful measurements: 

i) A logarithmic variation of the magnetization with field 
in the low- field two-dimensional vortex liquid below Tkt, 
with a coefficient set by the supcrfluid density and the 
temperature. Specifically, 

dM ^ tt Ps (T) f _ 2k B T \ 
d log if ~ # V nps{T))' 

Thus the theory predicts this simple relation between 
these directly measurable quantities, that can be checked 
in any sufficiently two-dimensional material where this 
regime should exist. 

ii) The double logarthmic variation of M in Eq. (JTJ at 
Tkt or at intermediate fields away from Tkt- Optimisti- 
cally, one might hope for a direct fit to this functional 
form. This field dependence also implies that in an M 
vs. log if plot, the curves for T < Tkt, though straight 
in the low-field regime, should exhibit an upward curva- 
ture at larger fields (see Fig. 0). 

iii) Correspondingly, plots of M(T) at different fields 
should exhibit a systematic drift of the "crossing point" . 
This is clear from plotting our results, as in Fig. 
which explicitly shows the downward creep of the "cross- 
ing point" towards Tkt from above as the field is de- 
creased within the two-dimensional vortex liquid regime. 
This trend definitely appears to be there in the recent 
BSSCO data of Wang, et al. 4 for underdoped and opti- 
mally doped samples. 

IV. SUMMARY AND FUTURE DIRECTIONS 

The mechanism of superconductivity in the cuprates 
remains one of the outstanding puzzles of the physics of 
correlated electrons. Nevertheless the proposition that 
aspects of their finite temperature behavior can be under- 
stood as consequences of sizeable thermal fluctuations of 
the superconducting order parameter has gained support 
in recent work. In this work we have examined the effects 
of such fluctuations in the two dimensional limit near the 
Kosterlitz-Thouless transition and presented an asymp- 
totically exact calculation of the magnetization in this 
vortex liquid state. The preliminary comparison with the 
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FIG. 3: Magnetization at 1, 2, 3, 4, 5, 10, 15, 20, 25, 30, 35, 
40, 45 and 50 Tesla as a function of temperature. Parameters 
are as in Fig. 1, in particular, Tkt = 80 K. 



highly anisotropic cuprates such as BSSCO is encourag- 
ing. 

On the theoretical side, it would be very useful to 
extend our calculation to higher fields, e.g., by keeping 
terms of higher order in vortex density in deriving flow 
Eqs. ©, © and (fTOft . The high-field behavior is one 
aspect of experimental magnetization data that does not 
appear to be captured well by our calculation. Strong sig- 
natures of superconducting fluctuations are also present 
in the Nernst coefficient. Indeed, Ref. |j has reported that 
the magnetization and the Nernst signal track each other. 
It would be desirable to have a theory of the Nernst ef- 
fect near the Kosterlitz-Thouless transition. We hope to 
report on this in the near future. 
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APPENDIX A: SUPERCONDUCTOR IN A FIELD 
TO COULOMB PLASMA 

For completeness, we sketch the mapping between 
a superconductor in a transverse field and a non- 
neutral Coulomb plasmaii. The Hamiltonian of a two- 
dimensional superconductor with a uniform magnetic in- 
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duction fioH c i <C B <C HoB C 2 can be approximated as 

2c 



(Al) 



where 9, A and = z • V x A are the usual phase, gauge 
and perpendicular magnetic fields. We also impose the 
constraint that the net vortex charge density matches 
B. Next, we explicitly separate the phase field into its 
longitudinal and transverse components by introducing 
a spin- wave field <f> and vortex charge density field n and 
Fourier-transform the Hamiltonian 



Hi 



sc 



q#0 



p s0 z x q 2 

= -5-|q<Pq + "q— 2~ I 

z (7 



-|q<^q 



PaO | ZXq, 2 



(A2) 

(A3) 
(A4) 



We then ignore the spin-wave part as it decouples from 
the rest of the problem. Transforming back into real 
space we arrive at the Hamiltonian of a non-neutral two 
component Coulomb gas 



H CG ^^f |d 2 xd 2 yMx)log 



(x-y) 2 



My),(A5) 



where Sn(y) = n(y) — B/(f>o. The range of integration 
above excludes the infinite self-interaction x = y of each 
vortex. One last remaining ingredient is the bare core 
energy E c q determined by the physics omitted in this 
derivation, it is usually taken as the energy cost of sup- 
pressing the order parameter inside the vortex core. So, 
finally, the Hamiltonian of a superconductor in the ex- 
ternal field is written Bsc — Bqg + E c qNt, which is a 
continuum version of Eq. (7). 
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